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Abstract 

The fluid thin film equation ht = — {h n h xxx ) x — a\ (h m h x ) x 
is known to conserve mass J hdx, and in the case of a\ < 0, 
to dissipate entropy J h?l 2 ~ n dx (see [8]) and the L 2 -norm of 
the gradient J h x dx (see [3])- For the special case of a\ = 
a new dissipated quantity J h a dx was recently discovered 
for positive classical solutions by Laugesen (see |15|). We ex- 
tend it in two ways. First, we prove that Laugesen's functional 
dissipates strong nonnegative generalized solutions. Second, we 
prove the full a-energy / (± h a h\ - {a+ ^_ 1 n h +i)™a+™-n+2) ) dx 
dissipation for strong nonnegative generalized solutions in the 
case of the unstable porous media perturbation a\ > and the 
critical exponent m = n + 2. 
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1 Introduction 

It is well known that analysis of the existence, uniqueness and regu- 
larity of weak solutions for nonlinear evolution equations relies heavily 
on a priori estimates. Often, the physical energy or entropy which 
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originate from the related model can provide non-increasing in time 
quantities. Unfortunately, it is far from obvious how to construct 
new non-increasing Lyapunov type functionals. A general algebraic 
approach to the construction of entropies in higher-order nonlinear 
PDEs can be found in [14J and can be applied to analyse thin film 
equations with stabilizing porus media type perturbations. In this pa- 
per, inspired by Laugesen's result [Tj)| on dissipation, we prove that 
the energy functional introduced in [T5J dissipates strong nonnegative 
generalized solutions. However, our method of the proof is only appli- 
cable to some subset of the Laugsen's dissipation region [15j (see the 
shaded area on Figure 1). 

We study the longwave-unstable generalized thin film equation 

h t = -{h n h xxx ) x - a x (h m h x ) x , (1.1) 

where h(x, t) gives the height of the evolving free-surface. The expo- 
nent n plays a stabilizing role due to fourth-order forward diffusion 
term and the exponent m plays a destabilizing role due to backward 
second-order diffusion term for the case when a\ > 0. This class of 
equations originates from many physical/industrial applications involv- 
ing air- fluid interface. For example: the case n = 1, m = 1 describes 
a thin jet in a Hele-Shaw cell |10| . the case n = 3, m = — 1 describes 
Van der Waals driven rupture of thin films [19], the case m = n = 3 
describes shape of fluid droplets hanging from a ceiling [llj, and the 
case n = 0, m = 1 describes solidification of a hyper-cooled melt 
(this is a modified Kuramoto-Sivashinsky equation) |3J. 

To prove that the nonnegativity property is preserved in nonlinear 
thin film equation h t = —(h n h xxx ) x for n > 1 (case a\ = 0) Bernis and 
Friedman [3] used set of dissipated and conserved quantities: mass 
conservation J hdx = M, surface energy dissipation ^ J h^dx < 0, 
and entropy dissipation -| J h 2 ~ n dx < 0. The new so-called /3-entropy 
J h 2 ~ n+l3 dx was introduced by Bertozzi and Pugh [5J and indepen- 
dently and simultaneously by Beretta, Bertsch, Dal Passo [lj to extend 
this result to n > 0. They also successfully used this new entropy to 
obtain exponential with respect to the L°°-norm convergence toward 
the mean value steady state solution. To analyse this convergence rate 
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in ii/^-norm for the special case n = 1, ai = Carlen and Ulusoy 
jl] used the dissipated energy j h a h\ dx constructed by Laugesen [TB] 
for classical positive solutions. Exponential asymptotic convergence 
toward the mean value was also studied by Tudorascu in [18]. This list 
of connections between new properties of solutions in thin film PDEs 
proved by means of newly discovered dissipated quantities is far from 
complete. 

In this paper we prove that there exists a subinterval J of —1 < 
a < 1 (I depends on n only) and a nonnegative strong generalized 
solution such that for any a E I the full a-energy 

\t) — J (| h a hi - ( a+m " 1 n +i)( a+m _ n+2 ) ) dx 
n 

dissipates. For the unstable porus media perturbation case ai > 
this dissipation is proven under the assumptions that the total mass 
of the solution is less than or equal to the critical one, m = n + 2 and 
domain Q is unbounded or h is compactly supported. For the stable 
case a>i < no such assumptions are needed. 

We proceed as follows. First, we show the dissipation for the clas- 
sical solutions of the regularized problem and then we take this dis- 
sipation to the limit. We prove dissipation of the full a-energy for 
positive classical solutions of the regularized problem for any value of 
the coefficient ai and without any additional assumptions about the 
total mass of the solution or its support. However our method of tak- 
ing the dissipation to the limit due to the Bernis-Friedman method of 
regularization requires additional conditions for the case a\ > 0. 



3 



2 Auxiliary results to generalized weak so- 
lutions 

We consider nonnegative weak solutions to the following initial-boundary 
problem: 

( h t + {h n h xxx + a x h m h x ) x = in Q T , (2.1) 

( P ) \ fK" a > *) = B(°» *) for t > 0, z = M, (2.2) 
[ h{0,x) = h o {x)^0, (2.3) 

where h = h(t,x), Q = (—a, a), Qt = (0, T) x Q, n > 0, m > 0, and 
a\ G R 1 . We define a generalized weak solution in the Bernis- Friedman 
sense (see, e.g. [HE]). 

Definition 2.1 (generalized weak solution). Let n > 0, m > 0, and 

a\ G R 1 . A generalized weak solution of problem (P) is a function h 
satisfying 

h G C^ 2 ' 1/8 (Q T ) n L°°(0, T; H\Q)), (2.4) 
h t eL 2 (0,T;(H\n))'), (2.5) 
h G Cti(9), h%(h xxx + ai h m - n h x ) G L 2 (T), (2.6) 

where 7 = Q T \ ({/i = 0} U {£ = 0}) and h satisfies \2. 1\) in the 
following sense: 

T 

(h t (-,t),cj ) )dt- jj h n (h xxx + a l h m - n h x )<f )x dxdt =0 (2.7) 

for all cf) G C 1 (Qr) with <p(—a, ■) = <p(a, ■); 

h(-,t) — > h(-,0) = ho pointwise & strongly in L 2 (Q) as t — )■ 0, (2.8) 

fc(-a, t) = fc(a, t) Vt G [0, T] and 0(-a, £) = g(a, t) (2.9) 
/or z = 1, 3 at all points of the lateral boundary where {h ^ 0}. 
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Because the second term of ( 12. 7p has an integral over CP rather 
than over Qt, the generalized weak solution is "weaker" than a stan- 
dard weak solution. Also note that the first term of (12. 7p uses h t G 
L 2 (0, T; (H 1 (Q))'); this is different from the definition of weak solution 
first introduced by Bernis and Friedman [3]; there, the first term was 
the integral of h<p t . The proof of the existence of generalized weak 
solutions follows the ideas of |3J HJ El El El EH- 

Let 

J if >-« / {-i. -•->). 



{(/3-n+2)(/?-n+l) 
zhiz-zttp-n=-l, (2-10) 
-biz if P-n = -2, 

(G^\z))" = z?- n , and G (z) := G °\z). 

Theorem 1. Let a\ G R 1 , n > 0; m ^ n/2 for a\ > 0, and m > for 
a-y < 0. 

(oj [Existence.] Let the nonnegative initial data ho G H 1 ^) satisfy 



J G {h {x))dx < oo, (2.11) 



and either 1) ho(—a) = ho{a) = or 2) ho(—a) = ho(a) ^ and 
|i£(_ a ) = |^(a) /ioWs fori = 1,2,3. T/ien /or some time T loc > 
there exists a nonnegative generalized weak solution, h, on Qt 1oc i> n the 
sense of the definition \2.1[ Furthermore, 

heL 2 (0,T loc ;H 2 (n)). (2.12) 

Let 

£ (T):= J {\hl(x,T) - ai D (h(x,T))}dx, (2.13) 



n 

where Dn(z) := -, 2 ™w" +2 — -ttt- Then the weak solution satisfies 

' (m—n+l)(m—n+2) J 



£ (T) + J J h n (h xxx + ai h m - n h x ) 2 dxdt ^ £ (0), (2.14) 



{h>0} 
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h n h 2 xxx dxdt^ const < oo. (2.15) 

{h>0} 

for all T < T\ oc . The time of existence, TJ oc; is determined by a\, 
/ ho, \\h 0x \\2, and J G (h ). Moreover, Ti oc = +00 for a\ < 0. 
(b) [Regularity.] If the initial data from (a) also satisfies 

G^\h (x)) dx < 00 



for some — 1/2 < /3 < 1, (3 7^ t/ien i/iere exists < < T; oc snc/i 
that the nonnegative generalized weak solution has the extra regularity 

h~ E L 2 (0, H 2 (n)) and h~ E L 2 (0, t£>; Wftfi)). (2.16) 

T/ie tzme of existence, , is determined bya±, J ho, ||/toai||2? an d 
J Gf\ha). Moreover, = +00 for a x < 0. 

There is nothing special about the time T\ oc in the Theorem [TJ In 
the case a\ > and n/2 < m < n + 2 (or m = n + 2 and M < M c ), 
given a countable collection of times in [0,7} oc ], one can construct a 
weak solution for which these bounds will hold at those times. Also, 
we note that the analogue of Theorem 4.2 in |3] also holds: there exists 
a nonnegative weak solution with the integral representation 

T 

J (h t (-,t), cj))dt + J J {nh n ~ x h x h xx (\) x + h n h xx <p xx ) dxdt (2.17) 

Q T 

- (H h m h x (j) x dxdt = 0. 



Qt 



3 Dissipation of energy for nonnegative weak 
solutions 

The main result of the present paper is the following 
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Theorem 2. Let a x G R 1 , 1/2 < n < 3; m ^ n/2 /or cii > 0, and 
m > for ai < 0, and 

£ { * ] {T):= J qh a h 2 x (x,T)-a 1 D (h(x,T))}dx, (3.1) 
h 

where Dniz) := 7 — ; ^ttt? " +2 r^r, Q^d £n = £n(T). Then there 

uv ' (»+m— n+l)(a+m— n+2) ' u \ / "\ ' 

exists a non-empty subinterval I (see /73]/ /or i/ie explicit form of the 
I) of < a < 1 for \ < n < 3, and 0/ § - n < a < for |<n<3 
snc/i that for any a G / £/ie nonnegative weak solution from Theorem^ 
satisfies the following estimates: 

(i) if a\ ^ inen 

4 Q) (TK4 Q) (0); (3.2) 

(ii) if a x > then 

E { a) (T) ^ £ { a) (0) + Ci J J h a+3m - 2n+2 dxdt for m> n + 2; (3.3) 

4 Q) (T) ^ £[, Q) (0) + T(C 1 M 2 " + " +4 + C 2 M a+3m - 2n+2 ) (3.4) 
/or m < n + 2 and a > 2n — 3m — 1; 

^ (T) ^ £^ a) (0) + C 3 T M a+3m ~ 2n+2 (3.5) 

/or m < n + 2 and 2n — 3m — 2 < a ^ 2n — 3m — 1; 

£[, a) (T) ^ £[, a) (0) + C 2 TM a+n+8 for m = n + 2 and < M < M c . 

(3.6) 

i/ere C2 = if Q is unbounded or h has compact support. 

Remark 3.1 (Extra Regularity). In particular, the extra regularity 
h~^~ G L°°(0, T; if 1 (f2)) follows directly from Theorem^ Hence, 
h^(-,T) G H 1 (Q) for almost all T G [0,T^ } ] and therefore (-,T) G 
C 1 / 2 ^) /or a/mosi a// T G [0,2^]. A ssume that Tq is chosen such 
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that h™? (-,To) G C 1 / 2 (fi) and h(xo,T ) = at some x G fi. Tnen 
iaere exists a constant L such that 

h^(x,T ) = \h^(x,T ) - h^(x ,T )\ < L\x - xol 1 / 2 . 



Hence h{x,T Q ) < L«+ 2 |x - x$\<*+ 2 , i. e. h(.,T) G C a + 2 (p) /or almost 



every T G [0, T, (/3)l 



Remark 3.2 (Rate of decrease). For a\ < and 1/2 <n< 3 we can 

generalize the results from {9, Theorem 1.1] in the following way: 

J h a h 2 x (x, t) dx < C{\ + for < a < 0, 
h 

whence \\h — < C(l + t)~2 /or any nonnegative strong solution 
h. Here C = C{a\,a, n, h, £q , (0)), and ft, = j^r||Ao||i- proof is 

similar to J9$. 

3.1 Regularized Problem 

Given 8, e > 0, a regularized parabolic problem, similar to that of 
Bernis and Friedman [3], is considered: 

( fk + (fse(h)(h xxx + ai D'^h)h x )) x = 0, (3.7) 

(Ps,e) I Q(-a, t) = g(a, t) for t > 0, % = M, (3.8) 

[ /i(a;,0) = /io,fc(a;), (3.9) 

where 

:= f £ (z) + 5 = + 6, D»(z) := (3.10) 

Vz G R 1 , e > 0, s ^ 4. The 5 > in (13~T0D makes the problem 
(13 .7ft regular (i.e. uniformly parabolic). The parameter e is an approx- 
imating parameter which has the effect of increasing the degeneracy 
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from f(h) ~ \h\ n to f £ (h) ~ h s . The nonnegative initial data, h , is 
approximated via 

h ,s £ e C 4+7 (ft), /jo,fc ^ ^o,<5 + e° for some < 6 < 

%^(-a) = %^(a)fori = M, ( 3 - U ) 
^o,fe - > strongly in as 5, e — > 0. 

The e term in (13. lip "lifts" the initial data so that it will be positive 
even if S = and the S is involved in smoothing the initial data from 
H\n) to C 4+ ^(fi). 

Sketch of Proof : By Eidelman [T2| Theorem 6.3, p. 302], the regu- 
larized problem has the unique classical solution hg E £ C^ 7 ' 1+7//4 (fi x 
[0, Taj) for some time t$ £ > 0. For any fixed values of 5 and e, by 
Eidelman [121 Theorem 9.3, p. 316] if one can prove a uniform in time 
an a priori bound \hs e (x,t)\ < A$ e < oo for some longer time interval 
[0, Ti OCt se] (Tiocjs > Tfe) and for all x G Q then Schauder-type interior 
estimates [121 Corollary 2, p. 213] imply that the solution hs E can be 
continued in time to be in C 4 ^ 7 ' 1+7//4 (f2 x [0, Ti OCi g £ ]). 

Although the solution h$ e is initially positive, there is no guarantee 
that it will remain nonnegative. The goal is to take 5 — > 0, e — > in 
such a way that 1) T; OCj< 5 £ — > T[ oc > 0, 2) the solutions /i^ converge 
to a (nonnegative) limit, h, which is a generalized weak solution, and 
3) h inherits certain a priori bounds. This is done by proving various 
a priori estimates for hs £ that are uniform in 5 and e and hold on 
a time interval [0,7] oc ] that is independent of 5 and e. As a result, 
{hse} will be a uniformly bounded and equicontinuous (in the C^/j* 
norm) family of functions in x [0,TJ o J. Taking 5 — > will result in 
a family of functions {h e } that are classical, positive, unique solutions 
to the regularized problem with 5 = 0. Taking e — > will then result 
in the desired generalized weak solution h. This last step is where 
the possibility of non-unique weak solutions arise; see [l] for simple 
examples of how such constructions applied to h t = — (\h\ n h xxx ) x can 
result in two different solutions arising from the same initial data. 
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3.2 Dissipation of energy for positive solutions 




Lemma 3.1. Let a belong to the full domain shown on Figure 1, and 

£^(T) := J {±h a h 2 x (x,T)- ai D £ (h(x,T))}dx, (3.12) 
n 

where D"(z) := z a D" £ {z). Then the unique positive classical solution 
h £ of the problem (Po, £ ) satisfies 

£( Q) (T) ^ £^(0) +/i J I h a ~ 2 f £ {h)D' £ \h)h A x dxdt+ 

Qt 

eh J h a - s ~ A f 2 £ {h)hldx + e 2 k 2 J h a ~ 2s - A f £ {h)hldx, (3.13) 
n n 

where h L = ki(a, n, s) are constants, and fi = (j,(a, ai) such that /i(0, a\) - 
and /i(a, ai) ^ for a± ^ 0. 
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Note that, although we use the same convenient notations intro- 
duced in [T5], the proof of Lemma 13.1 1 has essential differences from 
the proof of Theorem 1 of [15]. Indeed, we introduce new ideas in 
order to estimate the lower-order term in the equation (13. 7p . In par- 
ticular, the new quantity N is introduced, the quantity R is modified, 
and so are the terms involving the regularization parameter e in (I3.19p . 

Proof of Lemma \3. 11 To prove the bound (I3.13p . multiply (13. 7p with 
5 = by —^h a ~ 1 hl — h a h xx — aiD' £ (h), integrate over fl, use integration 
by parts, apply the periodic boundary conditions (13.81) . to find 



n 



f(a-l) 



Zi a) (t) = - J h a f £ {h){h xxx + ai D^h)h x ) 2 dx- 

j //" -hir (h)[h.,..,.. r + ai Di{h)h x ) dx- 

n 

2 a f h ^KK Mh)(hm + aim) K) dx . (3,4) 
n 

The equality (I3.14p can be rewritten as 

|£( Q )(t) = -R 2 -2aRS- f(a- 1) RL, (3.15) 
where the quantities 

R:= ({h a f £ {h))^ 2 {h xxx + ai D^h)h x )\ = \(h a f £ (h)y/ 2 (h xxx + ai D' £ \h)h x )) , 
S ■= ({h a - 2 f £ {h)) l ' 2 h x h xx \ = \{h a - 2 f e {h)f' 2 h x h xx ), 

l .= ({h^f E {h)Y' 2 hi\ = \(h^f £ (h)y/ 2 hi), 

N := ((h a - 2 f £ (h)D':(h))^ 2 h 2 x \ = \(h«- 2 f £ (h)D»(h))V 2 hl), 

each represent half of an inner product in L 2 (Q). We will need the 
following integration by parts formulas 

= -§(«- 3) J h a - 4 f £ (h)h 6 x dx-l J h a - 3 f' £ (h)h 6 x dx = 
n n 

- |(a + n - 3)L 2 - l e (s - n) J h ' 3 ' 4 f 2 (h)h 6 x dx, (3.16) 



n 
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RL = -(a+n-2)SL-e(s-n) J h a ~ s - 6 f 2 (h)h x h xx dx-3S 2 + ai N 2 = 

n 

\(a+n-2)(a+n-3)L 2 -3S 2 +a l N 2 -e{s-n) J h a ~ s ~ 3 f £ {h)h x h xx dx+ 

h 

\e{s-n){a+n-2) J h a ~ s ~ 4 f 2 (h)h 6 x dx = l(a+n-2)(a+n-3)L 2 -3S 2 + 
n 

a\N 2 + \e{s - n)(2a + 3n - s - 5) J h a ~ s - 4 f 2 (h)h 6 x dx+ 

h 

le 2 (s - nf J h a - 2s - 4 f £ (h)hl dx. (3.17) 



Here we use the auxiliary equality f £ (z) = nz 1 f E (z)+e(s—n)z ( s+1 ) f £ (z) 
Thus, from (I3.15P we have 

f t £, £ a \t) + fa(a-l)(s-n) 2 J h a - 2s - A f^h)h 6 x dx + fa(a-l)N 2 = 

a 

-R 2 -2aRS- f (a - l){a + n - 2) (a + n - 3)L 2 + f (a - l)S 2 + 

^a(a - l)(s - n)(s - 2a - 3n + 5) J h a ~ s ~ 4 f 2 (h)h 6 x dx. (3.18) 

h 

Our next step is to express ( 13. 18ft as the negative of a sum of squares 
to obtain the energy dissipation. To achieve this, we use ( 13. 16ft and 
(l3~TTj) to deduce that for all k <G R 1 , 

|£(«)(t) = -(R + a S + kL) 2 + 0(S + |(a+n-3) L) 2 + 7 L 2 + n N 2 + 
eh J h a - s - 4 f 2 (h)h 6 x dx + e 2 k 2 J h a - 2s - 4 f*(h)h 6 x dx, (3.19) 
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where 



kj_ = ^(s - n) I 5/c(a -s + 3n-5) + ^f{a - l)(a - 2a - 3n + 5) + 
(a + n - 3)(f (5a - 3) - 6«) ) , k 2 = \{s- n) 2 ( 4/c - a(a - 1) J , 



7 = k 2 - ^n(a + n - 3) (5(2 - n) + 3(a + n - 3)) - 

§(a + n - 3) (5(q - l)(ar + n - 2) - (5a - 3) (a + n - 3)) = 
k 2 - ^re(a + n - 3) (a - ^=±) - ^a(a + n - 3) (a - . 

Now we have to choose the parameter k in such a way that ^ and 
7^0. In this case, the parameter fi > for a\ > 0, and fi ^ for 
ai ^ 0. According to j 15 j . we can find k such that ^ 0, and 7^0 
when 1/2 < n < 3, see also Figure 1 where this region was computed 
numerically by Matlab (see |15j for the explicit form of the domain). 

3.3 Limit process in (13.131 ) 

Rewrite the integral JJ eh^~ s ~ 4 f 2 (h e )h^ x dxdt in the form 




Qt 




Qt 



Qt 



Using the Young's 



inequality 




(3.20) 



e 
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choosing p = — - > 1 and q = s n > 1 (=>- < a < s — n), we find 

Similarly, we deal with the integral e 2 f h £ ~ 2s ~ 4 f £ (h)h £X dx. Due to 

Q 

Lemma [A. II and (|2.15|) . ff h £ ~ A h^ x dxdt is uniformly bounded then 

Qt 

{kieh a - s -*f 2 {h e ) + k 2 e 2 h°- 2s - 4 f3{h £ ))h e EX dxdt ^ 

Qt 

Ce^ [f h n ~ A h% x dxdt ^Ce^, (3.21) 



where the positive constant C is independent of e. Letting e — > 0, from 
( 13.21 j) we obtain 

(heh«- s - 4 f 2 (h £ ) + k 2 e 2 h a £ - 2s - A f £ {he))hl x -> in L\Q T ) (3.22) 

for < a < s — n. 

Now, we show (13. 22ft for the case of a < 0. Rewrite the integral 
ff eh £ ~ s ~ A f 2 (h £ )h £X dxdt in the form 

Qt 

ehr s -*f s {h £ )h% x dxdt= J! ^^hr^dxdt. 

Qt Qt 

s— n , v i 

Using the inequality ( 13.201) with a = z p and b = 9 ' we obtain 

r z s - 2 < (p-l)f £ z°-i (p-l)f j(°-i)-ns- n ) 

C ( z s-n +£ \2 *5 „2 2(s-n) 2 p 2 fc ^ ' 

Z P ei 

choosing p = 2 f, s ~"^ 5 < 2 and o = tt, t^-t? — > 2 (=>- n > 

2 — a + /3), we find 



r ^ < ^-2 +a -^ 2/ i -2(n-l)-q+^ ^^ r ^. 
fc (z°-™+s) 2 ^ ^ 2(s-n) J I S-2+Q-/3 J fc ^ 
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where (3 G (—1/2,1) follows from ( I2.16p . Similarly, we deal with the 

integral e 2 J 7^~ 2s - 4 f*(h)h 6 EiX dx. Due to h G L°°(0, T; H 1 ^)) and 
n 

( I2.16p . ff h^~ 2 h\ x dxdt is uniformly bounded then 



It 



Ce " // hP- 2 hl x dxdt ^ Ce " " // h^- 2 hl x dxdt^ 



iT 

n — 2+a 



Ce — , (3.23) 

where the positive constant C is independent of e. Letting £ — )■ 0, 
we obtain f l3.22D for | — n < a < and | < n < 3. In view of the 
Lebesgue's theorem, we have 

h a - 2 j e {h £ )D" e {h £ )h\ x dxdt^ [ [ h a+m - 2 h x dxdt (3.24) 



Qt Qt 

if m > and a > -\ - m, due to h G L°°(0, T; and f l2~T6|) . 

Integrating ( I3.19P over the time interval, and letting e — > 0, in 
view of ( I3.22p and f )3.24p . we obtain ( 13.1 3p for some subinterval / for 
< a < 1 and | < n < 3 or for — 1 < a < and | < n < 3. 
Note that, the convergence on the left-hand side follows from Fatou's 
lemma and from the corresponding a priori estimate (see, for example, 

[31 EJECT]). □ 

3.4 Proof of Theorem [2] 

Taking the limit e — > we obtain 

4 a) (T)+7 / f h a+n - 4 h 6 x dxdt ^ £jj a) (0)+/z / / h a+m - 2 h x dxdt. (3.25) 



Qi 
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Now, we estimate J J h a+m 2 h x dxdt. Using the Holder inequality, we 

Qt 

obtain 



h a+m - 2 h 4 x dxdt < / ( / h a+n - A hidxY( ! h a+3m ~ 2n+2 dxY dt. 



q t on n 

(3.26) 

Applying Lemma IA.2I to v = h i with a = 6( - a " l ^ n n ^ +2 - > , d = 6, 
b = , 6 , „ < a (=>- a > 2n — 3m — 1), z = 0, and j = 1, we deduce 

//> q + 3m-2n+l p 3(2q + 5m-3n+4) 

h a + 3m-2n + 2 dx ^ j U ^ U hdxj ^ + 

n n n 

f f N«+3m-2n+2 3(2 Q + 5m -3n+4) / f , , - \ 

d 2 ( hdx) < ClAf f / h a+n ~ 4 hldxj + + 

a n 

C2M a+3m-2n+2 

Substituting (I3T27D in fl3T26|) . we find 

T 

h a+m ~ 2 h x dxdt <: Cl M 2 ° + " +4 / ( / h a+n -*h 6 x dx) a+n+7 dt+ 



o n 

T 



c 2 M a+3m ; 2 " +2 j ( | h a+n - 4 htdxfdt. (3.28) 



o n 

If m < n + 2 then, using Young's inequality, from (I3.28P we arrive 

at 

' r h a+m - 2 hidxdt 1 1 h a+n - 4 hldxdt+ 



C(e)r( Cl M 2Q " n+4 + c 2 M" +3m - 2 "+ 2 ). (3.29) 
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Substituting (13.291) in (I3.25p . and choosing e small enough, we obtain 

4 a) (T) ^ £i Q) (0) + T(C 1 M 2 " + " +4 + C 2 M a+3m ~ 2n+2 ) (3.30) 

for a > 2n — 3m — 1 and m < n + 2. Here C2 = if fHs unbounded or 
/i is compactly supported. In particular, if < a + 3m — 2n + 2 ^ 1, 
i. e. 2n — 3m — 2 < a ^ 2n — 3m — 1 then, using the Holder inequality 
and applying Young's inequality, from (I3.26P we obtain 

h »+m-2 h 4 dxdt ^ e J J h a+n ~ 4 hldxdt+ 
Qt 

C{e)\Vt\ 2n ~ 3m " 1 - a T M a+3m ~ 2n+2 (3.31) 

for 2n — 3m — 2 < a ^ 2n — 3m — 1. Substituting (I3.31[) in (I3.25p . and 
choosing e small enough, we obtain 

£ { a) (T) ^ £f, a) (0) + c 3 TM a+3m - 2n+2 . (3.32) 

If m = n + 2 then, using Young's inequality, from (I3.28P we deduce 



h a+m ~ 2 h x dxdt ^ Cl M 2 j j h a+n ~ 4 h x dxdt+ 

Qt 

e JJ h a+n ~ 4 hldxdt + C(e)T M a+n+s . (3.33) 

Qt 

Substituting (13.331) in (I3.25p . and choosing e enough small, we obtain 

£[, a) (T) < £^ a) (0) + C 2 T M a+n+8 (3.34) 

for a > -n-7,m = n + 2 and M ^ M c . Here C 2 = if Q is 
unbounded or /i is compactly supported. 
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Appendix A 

Lemma A.l. fl23 W Let Q C R N , N < 6, be a bounded convex 
domain with smooth boundary, and let n £ (2 — </l — ^g,3) /< 



for 



N > 1, and | < n < 3 /or N — 1. Then the following estimate holds 
for any positive functions v £ H 2 (Q) such that Vr • n — on <9fi and 
/ r n |VAr| 2 < oo: 



^ 6 K- 4 iv^r+^ 2 i^riv^r} ^ 

c{J V 6 v n \VAv\ 2 + J r n+2 |V^| 6 }, 

n W>o} 

where if £ C 2 (f2) zs an arbitrary nonnegative function such that the 
tangential component of V<£> egita/ io zero on <9f2 ; and tae constant 
c > is independent of v. 

Lemma A. 2. (116]) If Q C R a bounded domain with piecewise- 
smooth boundary, a > 1, b £ (0,a), d > 1, and ^ « < j, i,j £ N, 
taen iaere exist positive constants d\ and d 2 (d 2 = ifVtis unbounded) 
depending only on Q, d, j, b, and N such that the following inequality 
is valid for every v(x) £ W^ d {VL) fl L b {VL): 

l _|_ j l 

II^L«(Q) ^ rf l \\ D]V \\ 6 L^) W V W l L*U+ d 2 \H Lb{fl) , = \ N L _ 1 £ 

b N d 
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